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Recent advances in theory and computational experiments have shown the need to refine the 
previous categorisation of magnetic reconnection at three-dimensional null points - points at which 
the magnetic field vanishes. We propose here a division into three different types, depending on 
the nature of the flow near the spine and fan of the null. The spine is an isolated field line which 
approaches the null (or recedes from it), while the fan is a surface of field lines which recede from it 
(or approach it). 

So-called torsional spine reconnection occurs when field lines in the vicinity of the fan rotate, 
with current becoming concentrated along the spine, so that nearby field lines undergo rotational 
slippage. In torsional fan reconnection field lines near the spine rotate and create a current that is 
concentrated in the fan with a rotational fiux mismatch and rotational slippage. In both of these 
regimes, the spine and fan are perpendicular and there is no flux transfer across spine or fan. The 
third regime, called spine-fan reconnection, is the most common in practice and combines elements 
of the previous spine and fan models. In this case, in response to a generic shearing motion, the null 
point collapses to form a current sheet that is focused at the null itself, in a sheet that locally spans 
both the spine and fan. In this regime the spine and fan are no longer perpendicular and there is 
flux transfer across both of them. 



PACS numbers: 

I. INTRODUCTION 

Magnetic reconnection is a fundamental process of en- 
ergy release that lies at the core of many dynamic phe- 
nomena in the solar system such as solar flares, coronal 
heating events, geomagnetic substorms and flux trans- 
fer events. Reconnection in three dimensions has been 
shown to be completely different in many fundamental 
respects from the classically studied process in two di- 
mensions [1-3]. The main thrust of reconnection theory 
at present is to understand the different ways in which it 
may take place in three dimensions [e.g., the books 4, 5]. 
A key point is that in three dimensions reconnection oc- 
curs where a component of the electric field parallel to the 
magnetic field is present - and this can be in many differ- 
ent field configurations. For example, reconnection may 
occur either at null points [e.g., 6-13] or in the absence of 
null points at quasi-separatrix layers or hyperbolic flux 
tubes [14-23] or it may occur along separators that join 
one null point to another [7, 24-31]. 

Null points are common in the solar atmosphere [32- 
35] and are sometimes implicated in solar flares and coro- 
nal mass ejections [36-41]. Three-dimensional collapse of 
a null has been described [42-46] and stationary resistive 
flows near them have been modelled [47-49]. In partic- 
ular, for a linear null and uniform magnetic diffusivity, 
Titov and Hornig [49] discovered field-aligned flows when 
the spine current is small and spiral fleld-crossing flows 
which do not cross the spine or fan when the spine cur- 
rent exceeds a critical value. 

A three-dimensional null point possesses two different 
classes of field lines that connect to the null: for a so- 



called positive null point, a surface of field lines (called a 
fan by Priest and Titov [7]) recede from the null, while 
an isolated field line (called the spine of the null) ap- 
proaches it from two directions; for a negative null point, 
on the other hand the fan approaches the null, while the 
spine recedes from it. (For an alternative nomenclature 
see Ref. [50].) The different types of linear null were 
categorised by Parnell et al. [45]. The generic null in a 
potential magnetic field is an improper radial null, with 
the fan perpendicular to the spine and the field lines in 
the fan approaching or receding from essentially two di- 
rections (Fig. lb). A particular case is the proper radial 
null in which the field lines in the fan are radial (Fig. la). 
The effect of a current along the fan is to make the fan 
and spine no longer perpendicular (Fig. 3b), whereas a 
strong enough current along the spine makes the fan field 
lines spiral (Fig. 3a). 

There have been three steps towards categorising re- 
connection at a null point due to (i) analytical ideal mod- 
elling, (ii) kinematic resistive modelling and (iii) compu- 
tational experiments. The initial analytical ideal treat- 
ment by Priest and Titov [7] aimed to understand the 
types of ideal motions that arc possible in the environ- 
ment of a null point. They supposed that the nature of 
reconnection is determined to a large extent by the na- 
ture of the large-scale flows: they suggested that an ideal 
flow across the fan would drive spine reconnection, in 
which a current forms along the spine, whereas an ideal 
flow across the spine would drive fan reconnection with 
a strong current in the fan. They also proposed separa- 
tor reconnection with a strong current along a separator 
joining two nulls. 
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Since then, as we shall see in this paper, although be- 
haviour reminiscent of the early spine and fan models 
may be observed in certain limiting situations, recent 
numerical experiments have suggested different forms of 
spine and fan reconnection and also a hybrid spine-fan 
regime as being the generic modes that occur in prac- 
tice. However, the existence of separator reconnection 
has been well confirmed by a series of numerical exper- 
iments [25, 28, 31, 51] and its importance in the solar 
corona has been stressed [27, 30]. In addition, quasi- 
separatrix layer reconnection (called slip-running recon- 
nection by Aulanier et al. [40]) has been confirmed in 
numerical experiments [39, 52-57] and in bright point 
and flare simulations [17, 58-61]. 

Our aim here is simply to look more closely at the na- 
ture of reconnection at a 3D null point and to propose a 
new categorisation to replace spine reconnection and fan 
reconnection. In the next section it is necessary to sum- 
marise the main results from theory and computational 
experiments on null-point reconnection and to reinterpret 
them in the light of the new regimes of reconnection that 
we are proposing. In the following sections we consider 
in turn the properties of the three new types of recon- 
nection, namely, torsional spine reconnection, torsional 
fan reconnection and the most common regime spine-fan 
reconnection. 




II. THEORY AND NUMERICAL 
EXPERIMENTS 

A. Null Points 




FIG. 1: Field lines for (a) a proper radial null and (b) an 
improper radial null. 



The simplest linear null point (for which the magnetic 
field increases linearly from the null) has field compo- 
nents 



in Cartesian coordinates or 



So 



(1) 



iBR,B^,B,) = ^(i?,0,-2^) 

in cylindrical polars, so that V • B = identically, where 
J5o and Lq sue constant. The field lines are given by 



where c and k are constants. The z-axis is the spine and 
the a;y-plane is the fan. 

For this so-called proper radial null the fan field lines 
are straight (Figure la). It is a particular member (with 
a = 1) of a wider class of current-free improper radial 
null points (a 7^ 1) with curved fan field lines, having 
field components 

{B^,By,B^) = —[x,ay,-{a+ l)z\. 
^0 



This is the generic form for a current-free null, since the 
proper radial null is structurally unstable in the sense 
that it occurs only for a particular value of a, but for 
simplicity much of the theory so far has used a proper 
radial null. 

More generally, each of the three field components of 
a linear null may be written in terms of three constants, 
making nine in all. However, Parnell et al. [45] built on 
earlier work [62-64] and showed, by using V • B = 0, 
by normalising and by rotating the axes, that the nine 
constants may be reduced to four constants (a, 6, j±) 
such that 




1 kib-Jw) 




a 




where is the current parallel to the spine and j±/^J. 
is the current perpendicular to the spine. Furthermore, 
both nulls and separators are susceptible to collapse to 
form current sheets when the boundary conditions allow 
it [12, 24, 42]. 
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FIG. 2: Regimes envisaged from ideal motions: (a) Spine 
reconnection with a strong spine current driven by continuous 
motions across the fan. (b) Fan reconnection with a strong 
fan current and flipping of field lines above and below the fan 
produced by continuous motions across the spine. 



B. Kinematic Ideal Models 

The effects in the ideal region around a 3D null of 
steady reconnection were studied in the kinematic regime 
by Priest and Titov [7] extending earlier ideas [50] . They 
solved the equations 



and 



E + V X B = 



V X E = 



(2) 



(3) 



for V and E when B is given by Equation (1) and a variety 
of different boundary conditions are imposed. 

In particular, Eq. (3) implies that E = and then 
the component of Equation (2) perpendicular to B yields 



B • V$ = 0, 



(4) 



which, for certain imposed boundary conditions, may be 
integrated along field lines (characteristics) to determine 
the value of $ (and therefore E) throughout the volume. 
Then the component of Eq. (2) perpendicular to B de- 
termines the plasma velocity normal to B everywhere as 



V$ X B 

52 ■ 



(5) 



If a continuous flow is imposed across the fan (Figure 
2a), singularities in E and v are produced at the spine. 
Priest and Titov [7] speculated that this would produce 
a strong current at the spine in what they dubbed spine 
reconnection. They considered the effect of diffusion in a 
preliminary manner, but they were unable at the time to 



resolve the singularities at the spine. As an example, they 
considered flows with no (/)-component and an electric 
field of the form ~ VeBq sine/) giving rise to a velocity 



^EcfiL^z/Bo 

i?(i?2 +4z2)' 



E^L'^z/ Bq 



for which v±z is continuous at the fan z = 0, while v±r 
is singular at the spine R = 0. 

If, on the other hand, a continuous flow is imposed 
across the spine (Figure 2b), singularities are produced 
at the fan together with a strong flipping flow (that Priest 
and Forbes [65] had previously discovered). Priest and 
Titov [7] suggested that this would produce a strong cur- 
rent at the fan in what they dubbed fan reconnection. 
A particular example is given in terms of a; = x/L, 

y = y/L, z ~ z/i by a potential of the form $ = 
1 

VeBe[x'^z/ {A + y'^z)'^], which produces a flow field 

{v±x,V±y,V±z) 



(x^ -I- y2 + 4z2)(4 + y^zfl'^ 



2xyz{z^-l) 2{x'^ +Az^ + y'^z^) 



yl/2 



?l/2 



,{A + y'^z + x'^z)yz^' 



for which v±y is continuous on the planes z = ±1, while 
v±x and v±y are singular at the fan (z = 0). However, 
this analysis left open the questions as to whether it is 
possible to resolve the singularity and also whether these 
pure states are likely to be set up in practice. 



C. Kinematic Resistive Models 

The next step in the theory was to consider the effect 
in 3D of an isolated diffusion region where frozen-in flux 
breaks down and the induction equation is typically of 
the form 



dB 

'dt 



= V X (v X B) 



yyV^B. 



Reconnection in 3D is very different in many respects 
from that in 2D. 

In 2D, a diffcrentiable flux-transporting velocity w [66] 
satisfying 



dB 

'dt 



V X (w X B) 



always exists apart from at the X-point itself. This veloc- 
ity has a hyperbolic singularity at an X-typc null point, 
where the reconnection takes place. The magnetic flux 
moves at the velocity w and slips through the plasma, 
which itself moves at v. Furthermore, the mapping of 
the field lines in 2D is discontinuous at the separatrix 
field lines that thread the X-point. This mapping discon- 
tinuity is associated with the fact that field lines break 
and reconnect at one point, namely, the X-point. While 
they are in the diffusion region, field lines preserve their 
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connections everywhere, except at the X-point. Two flux 
tubes that move into the diffusion region break and rejoin 
perfectly to form two new flux tubes that move out. 

In 3D, surprisingly, none of the above properties carry 
over and so the nature of reconnection is profoundly dif- 
ferent [3]. First of all, a single flux tube velocity (w) 
does not generally exist [66, 67] since E ■ B ^ 0, but 
it may be replaced by a pair of flux velocities describ- 
ing separately what happens to field lines that enter or 
leave the diffusion region [3]. Secondly, the mapping of 
field lines is continuous if there is no 3D null point or 
separatrix surface. Thirdly, as they move through a 3D 
diffusion region, magnetic field lines continually change 
their connections. Fourthly, two tubes don't generally 
break and reform perfectly to give two other flux tubes: 
rather, when the two flux tubes are partly in the dif- 
fusion region and so are in the process of reconnecting, 
they split into four parts, each of which flips in a differ- 
ent manner, a manifestation of the continual change of 
connections. (Note that in general in 2D and 3D the flux 
velocity w is non-unique. We choose here to consider the 
case in which we select w by insisting that w = v in 
the ideal region. The crucial distinction is that in 2D a 
single w exists (and is singular), while in 3D reconnec- 
tion no single velocity w exists that satisfies (9) together 
with the constraint that w = v in the ideal region. See 
Refs. [66, 67] for further discussion.) 

The first attempt to model kinematically the effect of 
an isolated diffusion region was by Hornig and Priest [10] 
who set up a formalism and applied it to a case without 
null points. They solved 



E + vxB = 77j, 



(6) 



where VxE = 0, j = Vx B/^ and V • B = 0. The 
idea was to impose a sufficiently simple magnetic field 
that both the mapping and the inverse mapping of the 
field can be found analytically. Then, after writing E = 
V$, the integral of the component of (6) parallel to B 
determines $ everywhere as an integral 



$ = 



B 



along field lines, in terms of the values (<I>e) at one end of 
the field lines and the distance s along field lines. More 
simply in terms of a dimcnsionless stretched distance S 
such that ds/B ~ LodS/Bg, 



$ = 



Bo 



dS + ^e 



(7) 



One way of isolating the reconnection region in these 
kinematic solutions is by choosing a form of 77 that is 
localised. So-called pure solutions have <I>e = and pro- 
duce counter-rotating (or fiipping) flows of field lines that 
link the diffusion region. The rate of flux reconnection is 
calculated by evaluating the integral 



map 



dt 



(8) 





FIG. 3: The field near a null point with (a) uniform spine 
current and (b) uniform fan current. 



along a field line through the diffusion region [20, 68]. 
Then the flow normal to the field lines is determined by 
the component of Equation (6) perpendicular to B as 



(V$ - ?7 j) X B 
52 ■ 



(9) 



These solutions may be regarded as either kinematic (i.e., 
satisfying just the induction equation) or as fully dy- 
namic in the limit of uniform density and slow flow (since 
they also satisfy the equations V- v = and Vp = j x B). 

Pontin et al. [69] applied this formalism to determine 
the behaviour of the magnetic flux when an isolated dif- 
fusion region contains a spiral null point, i.e. a null with 
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current directed parallel to the spine line. The imposed 
magnetic field was 

[B^, By, B^) ^ — (x- ijoy, y + hjnx, -2z) 

or 

{Br, B^, B,) = ^ {R, \3oR, -2z) (10) 

in cylindrical polars, with the spine and cmrent both di- 
rected along the z-axis, where jo is a dimensionless cur- 
rent density. The diffusion region was assumed to be a 
cylinder of radius a and height 2b (Figure 3a). 

First of all, a pure elementary solution which describes 
the core of the reconnection process was obtained by set- 
ting the flow to zero outside the volume defined by the 
'envelope' (F) of flux that threads the diffusion region. 
Inside F the flow and flux velocities are purely rotational 
(i.e., in the ^-direction), so that there is no flow across ei- 
ther the spine or the fan. The reconnection rate is J E^^dl 
along the spine, and measures the rate of rotational mis- 
matching of the flux velocities of field lines entering and 
leaving the diffusion region. 

To this solution any ideal solution ($id) may be added 
and in particular they considered a stagnation-point flow 
of the form = (j)QXQyo, which brings flux into F and 
carries it out again. The result is a transition from O- 
type to X-type flow near the null when c/jq exceeds a crit- 
ical value. What this solution suggests, therefore, is that 
a type of spine reconnection with strong current along the 
spine direction is possible when there are twisting flows 
about the spine. This is quite different from the spine re- 
connection that was envisaged in Priest and Titov [7] and 
so here we propose to call it torsional spine reconnection 
and discuss its properties further in Section III. 

Next. Pontin et al. [70] applied the same approach to a 
diffusion region (D) containing a null point having a uni- 
form fan- aligned current (Bojo/il^Lo)) in the x-dircction 
and field components 

{B.j„By,B^) = —{x,y - joz,-2z). 

^0 

The diffusion region was assumed to have the shape of a 
disc of radius a and height 26 (see Fig. 3(b)), inside which 
the magnetic diffusivity decreases smoothly and mono- 
tonically from the null to zero at its boundary. Outside 
D it vanishes. 

The resulting plasma flow was surprisingly found to 
be quite different from the fan reconnection of Priest and 
Titov [7] , since it is found to cross both the spine and fan 
of the null. Field lines traced from footpoints anchored 
in the fan-crossing flow are found to flip up and down 
the spine, whereas those that are traced from the top 
and bottom of the domain flip around the spine in the 
fan plane, as envisaged by Priest and Titov [7]. The 
reconnection rate is again given by an integral of the 
form (8), this time along the fan field line parallel to the 



direction of current flow (here the x-axis). For such a 
mode of reconnection this expression can be shown to 
coincide with the rate of flux transport across the fan 
(separatrix) surface [70]. 

It is possible to find a solution that has similar field 
line behaviour to the pure fan reconnection envisaged 
by Priest and Titov [7], with flow across the spine 
but not the fan, by adopting instead a field of the 
form (Bq/ Lo){x,y — Jqz^ / Lq, ~2z) with a fan x-current 
3Bojoz^ /{fiLg) (see Ref. [70]). It is also possible to model 
pure spine reconnection with flow across the fan but 
not the spine by considering [Bq/ Lq){x, y, joy^ /Lq — 2z). 
with a fan x-current SBojoV^ /{ij-Lq). Both of these fields 
have a vanishing current at the null. However, a key 
property of a null point is the hyperbolic field structure, 
which tends to focus disturbances and thus generate non- 
zero currents at the null for the primary reconnection 
modes. The above pure spine and fan solutions should 
therefore not be considered as fundamental or primary 
reconnection modes but as secondary reconnection modes 
in the sense that the current vanishes at the null. 

It has been suggested that solutions for spine reconnec- 
tion in incompressible plasmas [47] may not be dynam- 
ically accessible, and while incompressible fan solutions 
[6] are dynamically accessible [8, 71-73], this breaks down 
when the incomprcssibility assumption is relaxed [71]. It 
turns out that the generic null point reconnection mode 
that is observed in numerical experiments in response to 
shearing motions is one in which there is a strong fan 
current with flow across both spine and fan, and which 
is in some sense a combination of the spine and fan re- 
connection of Priest and Titov [7] . We propose here to 
call it spine-fan reconnection and discuss its properties 
further in Section V. 



D. Numerical Experiments 

Several numerical experiments have been conducted in 
order to go beyond the constraints of analytical theory 
and to shed more light on the nature of reconnection 
at a 3D null. The aim was also to see whether the 
types of reconnection envisaged qualitatively could in- 
deed take place in practice and to discover whether any 
other regimes are possible. 

First of all, Galsgaard et al. [74] investigated propa- 
gation of a helical Alfven wave towards the fan plane, 
launched by a rotational driving of the field lines around 
the spine. This led to the concentration of current in the 
fan plane and suggests the possibility of torsional fan re- 
connection which we shall propose in Section IV. (For 
highly impulsive driving coupling to a fast mode wave 
that wraps around the null was also observed.) On the 
other hand Pontin and Galsgaard [75] used a resistive 
MHD code to show how rotational disturbances of field 
lines in the vicinity of the fan plane can also produce the 
strong currents along the spine that are symptomatic of 
torsional spine reconnection (Section III). 
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FIG. 4: (Color) (a) A shearing motion of a spine that is situ- 
ated on the z-axis. (b) The resulting collapse of spine and fan 
to form spine-fan reconnection, showing the current-density 
contours (colour) and flow velocity (white) in the x — plane 
[After 77]. 



Then Pontin and Craig [12] used an ideal Lagrangian 
relaxation code to follow the formation of current sheets 
by the collapse of a line-tied 3D null in a compressible 
plasma. This was a result of the focussing of externally 
generated large-scale stresses in the field in response to 
an initial shearing of either the spine axis or fan plane. 
Building on a previous linear theory by Rickard and Titov 
[76], they found that locally the fan and spine collapse 
towards each other to form a current sheet singularity. 
This was followed up by Pontin et al. [77], who used a 
resistive MHD code to investigate the formation and dis- 
sipation of the current sheet in response to shearing of 
the spine, as shown in Figure 4. The results support 
the idea of spine-fan reconnection in which current con- 
centrates around the null (in a sheet spanning the spine 
and fan). Including compressibility does not affect the 
results qualitatively, except that in the incompressible 
limit the spine-fan current is found to reduce purely to a 



fan current [71] with behaviour closely resembling earlier 
fan reconnection models [6, 7]. So pure fan reconnection 
can be either an incompressible limit of spine-fan recon- 
nection or, as we have seen in Section II C, the result of 
a secondary fan current which vanishes at the null. 



III. TORSIONAL SPINE RECONNECTION 

The type of reconnection set up at a 3D null depends 
crucially on the nature of the flows and boundary con- 
ditions that are responsible for the reconnection. Let us 
suppose first that a rotation of the fan plane drives a 
current along the spine and gives rise to torsional spine 
reconnection, as sketched in Figure 5a. The nature of 
the reconnection is that in the core of the spine current 
tube there is rotational slippage, with the field lines be- 
coming disconnected and rotating around the spine (see 
Pontin and Galsgaard [75]): Figure 5b shows on the left 
side a particular magnetic field line and its plasma ele- 
ments at t = to; in the upper part of the figure (above 
the shaded diffusion region) this field line and its attached 
plasma elements rotate about the spine through positions 
at times ti, t2 and t^; in the lower part of the figure (be- 
low the diffusion region) the plasma elements that were 
on the field line at to rotate to positions at ti, t2 and 

that are on different field lines. A steady kinematic 
solution may be found following the approach of Sec- 
tion lie. The electric field may be written as the sum 
(E = V$ = V^nid + V<l>jd) of a nonideal pure (elemen- 
tary) solution satisfying 

V$™d + v„,d X B = r/V X B, 

and an ideal solution satisfying 

V4>,d + VidxB = 0. 

Consider a spiral null point (Equation 10) and suppose 
the diffusion region is a cylinder of radius a and height 
2b and that the magnetic diffusivity has the form 77 = 
?7o/(i?, z), where /(O, 0) = 1 and f{R, z) vanishes on the 
boundary of the diffusion region and outside it. 

The field lines for this spiral null may be obtained by 
solving 



dR L^Br d(p ^- dz 

dS^-B^'^' ^dS~-'^'^' dS 



-2z. 



Suppose we start a field line at the point {R, (f), z) 
(i?o, 4>Q,b) a,t S = Q. Then the field line equations are 



R^ Ra e 



s 



= 6 



+ \ JO S. (11) 



These give a mapping from an initial point (i?Oi 0Oi b) to 
any other point {R, (p, z) along a field line. The inverse 
mapping is 

-,-S A — A 



where S 



Rn = Re" 



ilog(z/&). 



jo S"- 



(12) 
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FIG. 5: (a) A rotational motion of tiie fan (open arrows) driv- 
ing torsional spine reconnection with a strong current (solid 
arrows) along the spine, (b) Rotational slippage of fields en- 
tering through the top of the diffusion region on a curved flux 
surface, showing as solid curves the locations of the plasma 
elements aX t = ti, t = t2, t = t^, that initially {t = to) lay on 
one field line, (c) The reconnection rate measures a rotational 
mismatching of flux threading the diffusion region, namely the 
difference between the rates of flux transport through surfaces 
A and B. 



A. Pure Non-Ideal Solution 

The pure elementary solution describes the core of the 
reconnection process. It is obtained following Refs. [10, 
69] by solving E + vxB = 77 j, with V x E = 0, 
j = V X B//1 and V • B = 0. Thus we write E = V^nid 
with ^nid given by Equation (7) and set <I>e = so that 
the flow vanishes outside the diffusion region. Inside the 
diffusion region the flow and flux velocities have no com- 
ponent across either the spine or the fan. For the spiral 
magnetic field {Bn,B^,B^) = {Bq/ Lo){R, ^joR, -2z) 
and the mapping (11), ^nid becomes 

^,iid = -^nido y 77/770 e'^-^dS", 

where 4>„ido = 2i?ofejo77o/(Mio)- Then, once a form 
for T] is assumed, this may be integrated to give 
^mdC'S', i?0: 0o)- After using the inverse mapping (12), 
we can then deduce ^nid{R,4>,z) and therefore E and 
v_L everywhere. 

If a diffusion region is isolated, a change of connectiv- 
ity of field lines may be studied, by following field lines 
anchored in the ideal region on either side of the diffu- 
sion region. A diffusion region is in general isolated if ?7j 
is localised in space. In practical cases in astrophysics, 
this is likely to be mainly because j is localised but in 
addition sometimes because as a consequence rj is also 
localised. Some numerical simulations have a localised 77, 
whereas others have a uniform rj or a purely numerical 
dissipation. But the important feature in all these cases 
is that the product rjj is localised. Now, in each of our 
solutions below, we follow Refs. [10, 69, 70] in choosing 
a spatially localised 77j by imposing a spatially localised 
resistivity profile together with a j that is not localised. 
The reason for doing this is to render the mathematical 
equations tractable, since we have not yet discovered a 
way to do so with a localised j. The quantitative spa- 
tial profiles of physical quantities will depend on the 77 
profile, but the qualitative topological properties of the 
field line behaviour in such models are expected to be 
generic and independent of the particular profile chosen 
for 77. Indeed, the topological properties of the reconnec- 
tion models of Refs. [10, 69, 70] have been verified by the 
numerical simulations [53, 75, 77]. 

There are four regions with different forms for ^nid, as 
illustrated in Figure 6, which shows a vertical cut in the 
first quadrant of the i?z-plane. In region (1) threaded by 
field lines that enter the diffusion region (shaded) from 
above, we assume ^nid{R,z) = 0, so that there is no 
electric field or flow. The same is true in region (2) which 
lies above the flux surface zR^ = ba^ that touches the 
upper corner (a, b) of the diffusion region. We calculate 
below the forms of ^nid{R,z) in the diffusion region (3) 
and in the region (4) threaded by field lines that leave 
the diffusion region through its sides. 

For example, let us assume that rj vanishes outside the 
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FIG. 6: The projection of magnetic field lines and the diffu- 
sion region in the first quadrant of the R-z plane, showing 4 
different regions (f)-(4) in which ^nid{R,z) is calculated. A 
magnetic field line whose projection intersects the top of the 
diffusion region in T{R,b) and the side in Q{a,Z3) contains 
typical points P{R, z) inside and beyond the diffusion region. 
The bounding field line zi?^ = hc^ is shown dashed. 



diffusion region (D) and that inside D it has the form 

-2 " 



^7 = 1 - -T 



z 



which peaks at the origin and vanishes on the boundary 
of D. First, we use the mapping (11) to substitute for 
R and z, and integrate with respect to S from the point 
T{R, b) on the top of D to the point P{R, z) inside D 
(Figure 6). Then wc use the inverse mapping (12) to 
replace i?o s-nd S*, and finally we obtain the potential 
throughout D (region (3) in Figure 6) as 



idO 



1 - 



- 1 



z 
b 

a4 



R^ 

a* 

2 



Z 
&3 



z z 

b^¥ 

3\1 



This then determines the components of the electric field 
(E = 'V^nid) everywhere in D as 



2^md0R^ 



dR 



2z2 

62 



id 



dz 



nidO 

~2b~ 



Rr__^ 



2 4zi?4 3z^R-i 



ba^ 



b^c 



In order to find ^nid{R, z) in region (4) of Figure 6, 
we start with the values of ^nid at the point Q(a, z^) 
on the side of the diffusion region (Figure 6) and then 
calculate ^nid at any point P{R, z) that lies on the same 
field line in region (4) to the right of Q. Thus, after 
putting {R, z)~{a, Zs) in the expression (13) for $ that 
holds in the diffusion region, we obtain 



^nid{a,Zs) = f{Zs) 



nidO 



62 53 



31 



(14) 

Since ideal MHD holds in region (4), ^nid{R,z) is con- 
stant along the field line {zR^ = ZgO^) joining Q to P, 
and so the value of ^nid at P is simply 

/zi?2 



nidO 



1 _ Zi?2 

3 ~ b~^ 



z^lV_ 



R^ 



63 



.(15) 



The solution for z < can be obtained in a similar man- 
ner by integrating from z = —6. 

We may now make various deductions from the solu- 
tion. The reconnection rate depends on the form of rj 
and is given in order of magnitude by 



Ei\ ds 



2Eo 6, 



where Eq is the electric field at the centre of the diffusion 
region and 26 is the dimension of the diffusion region 
along the magnetic field direction. In our example, Eq = 
£:^(0,0,0) = $™do/(26) = Tjjo, where jo = joBo/(m^o) 
is the value of the current at the origin, and along the 
spine (13) implies that 



£^.(0,0, z) 



nidO 



26 



1 



and so the reconnection rate becomes more accurately 

fb 



E;,{0,0,z) dz = ^ Eob 



-b 



(16) 



The other feature that wc can deduce from the elec- 
tric field components is the perpendicular plasma veloc- 
ity given by Eq. (9). In particular, on the fan plane 
(z = 0) inside D, En = 0, E, = ($™do/26)(l + RVa^), 



VJz = ($„»do/26)(l - Rya'^) and Br = BqR/Lo so that 
(13) there is a rotational component given by 



V4, 



(E, - mz)Bn 
52 



i?3 



where vq = ^nidoL^ / \baBo{l + \ Jq)]- The nature of the 
fiow becomes clear if we subtract a component parallel 
to B in order that Wz = (we are free to do this since the 
component of v parallel to B is arbitrary in the model). 
After doing this we find that vr vanishes, leaving v ~ 
(0, 0), i.e., the flow corresponds to a pure rotation (as 
in the solutions of Refs. [10, 69]). 
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B. Extra Ideal Solution 

To the above pure diffusive solution any ideal solution 
may be added satisfying E + v x B = and V x E = 0, 
for which the potential (<f>id) satisfies 

B ■ V$,d = 0. 

Thus, once the functional form ^id{RQ,(j)Q) is chosen at 
the points {Ro,(j)o,b) on z = b, say, that form of 
is constant along field lines given by the mapping (11). 
The resulting variation of <i>jrf(i?, (/i, z) throughout space 
is given by substituting for Rq and (f>o from the inverse 
mapping (12). 

As an example, suppose 



Rr, 




on the plane z = b. Then throughout the volume we find 



^^d{R,(t>:Z)^<^^M 2 



R^z 
a^b' 



which implies electric field components 



do 



2Rz, 



idO 



Then the plasma velocity components follow from 
E X B/B^ as 



where a 



1 



^ ^^doLo i-fjoR\R' +ARz^j„R^z) 
a^bBo {a^R^ + 4^2) 

(17) 

jJq. In particular, we notice that the 
flow vanishes on the spine i? = 0, and that in the fan 
z — there is a rotational flow that linearly increases 
with distance Vci,{R,(f>,0) = —^idoLoR/{a'^bBoa'^). 

The reconnection of field lines takes the form of a rota- 
tional slippage. Field lines entering the diffusion region 
have a flux velocity Wi„ = — V'S/n x H/B^, while those 
that leave it have a flux velocity Wout = — V$out x B/B^. 
$i„ is obtained by integrating along field lines that enter 
from the ideal region on one side, while ^out is obtained 
by integrating backwards along field lines that leave from 
the other side. The rate of slippage between inward and 
outward flux bundles is given by Aw = Wout — ^in and 
represents the rate of reconnection, which we have evalu- 
ated directly above in Equation (16). This reconnection 
rate, obtained by integrating along the spine, mea- 
sures the difference between the rates of ffux transport 
across surface A and surface B in Fig. 5(c). 

Note that the extra ideal solution does not change the 
rate of relative slippage. However, it does allow for differ- 
ent external conditions, such as rotation above and below 
the diffusion region in the same or opposite senses. To 
see the effect of a non-rotational ideal flow see Ref. [69] . 
In the solution given above the physical quantities E and 
V are continuous but not diffcrentiable at the boundary 



FIG. 7: A rotational motion of the spine (open arrows) driving 
torsional fan reconnection with a strong current in the fan and 
slippage of field lines (solid arrow). 



between regions (3) and (4). This is a sacrifice made for 
tractability and pedagogic purposes. For a solution with 
diffcrentiable physical quantities, see Pontin et al. [69]. 

In the above, the diffusion region was imposed to be a 
cylinder whose width (a) and height (26) are parameters 
of the solution. The formation, in a self-consistent fash- 
ion, of such a cylindrical diffusion region was observed 
in the simulations described by Pontin and Galsgaard 
[75]. In one of their simulations they imposed a twist- 
ing perturbation of the magnetic field in the vicinity of 
the fan plane. As the disturbance propagated inwards 
towards the null, it was dominated by a helical Alfvcnic 
wave - travelling along the field lines and thus stretching 
out along the spine. The result was a tube of current 
focussed around the spine, giving a large aspect ratio to 
the diffusion region (6 3> a). During the process of tor- 
sional spine reconnection the narrowing and elongation 
of the current tube is likely to continue until the rota- 
tional advection that twists the field and intensifies the 
current is balanced by the rotational slippage. 



IV. TORSIONAL FAN RECONNECTION 

Now suppose that we rotate the field lines near the 
spine in opposite directions above and below the fan. 
Then a current builds up in the fan. Within the fan 
current sheet, field lines experience rotational slippage 
[74, 75], in the opposite sense above and below the fan, 
in what we propose to term torsional fan reconnection 
(Figure 7). Again there is no fiow across either spine or 
fan. 

The counter-rotation (above and below the fan) of the 
region around the spine builds up a double-spiral struc- 
ture near the null point, with a current that possesses 
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two components: an axial component that reverses sign 
at the fan plane and a radial component. A counter- 
rotating part to the diffusion velocity {rjjuBz) is set up 
in the (;!i-direction that reverses sign at the fan. 

In order to model such reconnection, we consider what 
we term a double- spiral null point with field components 



Bn I - z 

^0 



2M+lj^N-l 



— , -2z (18) 



where M and N are positive integers and the correspond- 
ing current components are 



{-{2M + l)z^''R^-\ 0, 7Vz2M+i^w-2) ^ 

An alternative solution to the one presented below is out- 
lined in the Appendix. 

The field line equations for a mapping from an initial 
point (i?o, 00, b) to any other point (i?, 0, z) are 



i? = i?o e^, z^b e"2^, 



N-2 



and the inverse mapping is 
Ro = R e^^, 



1 - e 



-{4M-N+4}S 



2jo Ro 



N-2 



4A/ - iV + 4 5^-2 
-ilog(z/6). 



,-(4M-Ar+4)5 



where S 

Let us follow the approach of Section IIIA and calcu- 
late the pure non-ideal solution. We shall assume the 
diffusion region to be a disc of radius a and height 26, 
with the same diagram as before (Figure 6), except that 
the diffusion region is now expected to be in the shape 
of a thin disc (with b a) rather than a thin tube (with 
b ^ a). Assuming, as before, that ^{R,z) vanishes in 
regions (1) and (2), we evaluate it in region (3) by in- 
tegrating from a point T{R, b) on the top of the disc to 
a point P(i?, z) inside the diffusion region. After using 
Equation (7) and the mapping and setting <I>e = 0, the 
expression for the potential at P(i?, z) then becomes 



$ = -$™.o/^((2M-fl):^e 



'('iM-N)S 



+2iv|2_e-(4M-^+6)s ^s. (19) 

We adopt the following general form for the magnetic 
diffusivity inside the diffusion region (D) 



6" 



which peaks at the null point and vanishes on the bound- 
ary of D when to and n are positive and n is even. After 



substituting into (19) and using the mapping and inverse 
mapping, we find the potential throughout the diffusion 
region. In particular, it transpires that an important con- 
straint on the constants M, iV, to and n is that Ez be 
finite and continuous at the fan plane. As an example, 
one set of such constants that works is A/ = 2, iV = 6, 
TO = 4 and n ~ 2, for which 



<^md{R,z) = 



■ —^nidO 1 


z2i?4 5 


z^R^ 


bz-^R^ 
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369 


2610 


5z^^R^^ 


5z6i?i" 


z8i?4 
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26i2a4 
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3z^R^ 


3z^R^ 1 
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268a-i 
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26i2a4 J 



(20) 



The corresponding components of electric field are 



Er = 



nid 



nidO 



OR 



2z2i?3 10z3i?5 



Ibz^R^ bQz^R^ 5z^R^ 25z6i?9 
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~W ^ ~68 



69 



4 Rio 



'^z^R^ Ibz^R 
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2hz-^R 

b^^a^ ^ "611 511^ 
8z^i?4 10z3i?i" 9z5i?4 



6" a469 69 

&z^R^ 18z^R^ l2z'^R^ 



b'^a^ 



In order to find $(i?, z) in region (4) of Figure 6, as 
before, we calculate its value at any point Q{a, Zg) where 
a field line leaves the diffusion region, and then project 
that value along that field line. Thus, after putting 
(i?, z) = (a, Zs) into Equation (20), we obtain 

$md(a, 2s) = f{Zs) = -^nidO { ^^^ItTT + 



266 



369 



2610 611 6^12 2612 512 
5z!a6 3z6a4 3zW^ Sz^a'^ 3z^a^ 



2610 



261 



268 



610 



261 



(21) 



Since ideal MHD holds in region (4), ^nid{R,z) is con- 
stant along the field line (zi?2 = Zgu'^) joining Q to P, 
and so the value of ^nid at P is simply 



<^,m{R,z) = / 



Zi?2 
2~ 



nidO 



z'^R'^ hz^R^ 



a- / 1^ 266 359 

5z4i?8 5z5i?10 5z6i?12 5z6i?12 ^8i?16 



6a66i2 2a66i2 ai26i2 

5z4i?8 3^6^12 3^4^8 3^6^12 3^8^16^ 



2a26io 2a86io 20*68 ^s^io 2ai26i2 



(22) 



11 



The electric field components vanish in both the spine 
and the fan but are strong just above and below the fan, 
which is where the reconnection of field lines occurs by 
rotational slippage in a similar fashion to torsional spine 
reconnection. Near the spine and fan we have to lowest 
order in R and z 



The reconnection rate is the maximum value of 
(/ E\\ ds) along any field line {B?z = i?Q6), each of which 
enters the diffusion region from above at T{Ro,b) and 
leaves at Q{a,Zs). Along such field lines the integral is a 
function of Ro/a and b/a, namely: 



E\\ ds 
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E B 
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12z^R^ 
l2z'^R^ 
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E B 
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5z^R^ bz^R^ 
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12z6i?3 



69 
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nidO 



R^ 
264 



366 268 ^4 
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When i?o a >• 6 for a slender disc-shaped diffusion 
region, this reduces to 



E\i ds = — 



264 



9Rl 9i?J 

^8 " 



12 



12 



nl6 
„16 



If a and 6 are held fixed and Rq is varied, the maximum 
value of this occurs at Rq 0.90a, giving a reconnection 
rate of 



Eu ds 



0.9 ^ 



nidO 



64' 



(23) 



As for torsional spine reconnection, the reconnec- 
tion rate is proportional to the potential ^nido = 
2i3o6jo'7o/(Mio)i but in this case, as well as being pro- 
portional to the current density jo and diffusion region 
height (6), it also depends on its aspect ratio (a/b). 

Again, as before, a wide range of ideal solutions {^id) 
may be added to the diffusive solution. Thus, if for exam- 
ple, ^id{Ro,(j)o) = *ido-Ro/a" on the top (z = 6) of the 
diffusion region, the fact that it remains constant along 
field lines {R^z = R^b) determines 



from which the electric field components can be deduced. 



For instance, the case n = 4 gives an electric field of 

iER,E^,E^) = ^{2R'z^0,R^z), 

which implies a plasma velocity with components normal 
to the magnetic field of 

{vR, v^, w^) = {-E^B^, E^Br - ErB,_, ErB^) 

= {-^~J,z'R\ i2R'z + 8R'z')b^8joz'R') , 

g{R,z) 

where giR,z) = [R"^ + AR^° z^°j^/{b^^) + 4z^]Boa%^^ / Lq. 
In particular, it gives a rotational component (w^) that 
is odd in z and so represents the kind of counter-rotation 
that is typical of torsional fan reconnection. 



V. SPINE-FAN RECONNECTION 

In general, if the driving motions tend to shear a null 
point rather than rotate it, then the result will be spine- 
fan reconnection. A shear disturbance of either the spine 
or fan of the null will tend to make the null 'collapse'. 
That is, the resulting Lorentz force acts to increase the 
displacement, just as at a 2D null [see Ref. 12] and as 
at a separator [25, 26]. This collapse is opposed by the 
line- tying at the boundaries, and what occurs is that the 
shear distortion tends to focus in the weak field region in 
the vicinity of the null point, forming a localised current 
sheet [12, 77]. 

What distinguishes spine-fan reconnection from the 
other null point reconnection modes is that flux is trans- 
ferred across both the spine and fan. Furthermore, the 
current concentration is in the form of a localised sheet 
that is inclined at an intermediate angle between the 
spine and fan - indeed the current sheet contains part 
of both the spine and the fan (see Figures 4b and 8a). 
As mentioned above, the reconnection rate for this mode 
of reconnection is obtained by integrating E^\ along the 
fan field line with the maximum value of J E^^ds. By the 
symmetry of the simple models described herein, this is 
the field line parallel to the current orientation at the 
null (perpendicular to the applied shear). The reconnec- 
tion rate thus obtained measures exactly the rate of flux 
transport in the ideal region across the fan separatrix 
surface. To illustrate the properties of this mode of null 
point reconnection, we describe here briefly the results of 
resistive MHD simulation runs (see Pontin et al. [77] for 
an initial description). 

In the simulations, a shear velocity is prescribed on the 
(line-tied) z-boundaries, which advects the spine foot- 
points, see Fig. 4(a) (the results are qualitatively the 
same if the fan is distorted instead). The current sheet 
that forms in response to the shearing is localised in all 
three directions about the null. However, in the plane of 
the applied shear (perpendicular to the current orienta- 
tion at the null) the magnetic field and current patterns 
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FIG. 8: (a) The magnetic structure of the field in spine-fan 
reconnection, showing the field lines and the (shaded) diffu- 
sion region, (b) The corresponding motion of flux across both 
the spine and fan (large light arrows). The current sheet is 
shaded (with the part below the fan having a lighter shading 
than the part above) and contains a current flowing in the 
a;-direction (large dark arrows): its width is I, its total length 
Ltot in the yz-plane, and its length Lc common to spine and 
fan. 



have a similar appearance to a 2D X-point configura- 
tion. As one moves away from the null in the fan along 
the direction of current flow, the magnetic field, strength 
parallel to the current (sometimes known as the 'guide 
field') strengthens, while the current intensity weakens ~ 
see Fig. 8. 

The boundary shearing velocity is ramped up to a con- 
stant value (vq) at which it is held until being ramped 
down to zero, at t = t = 3.6 (space and time units in the 
code are such that an Alfven wave would travel one space 
unit in one time unit for uniform density and magnetic 
field |B| = l,yO = 1). The resistivity is uniform. Current 
focusses at the null during the period when the driving 
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FIG. 9: Scaling with the driving velocity vo of (top left- 
bottom right) Ltot, Lc, I and (see Figure 8b for notation). 



occurs, and when the driving ceases both the current and 
reconnection rate peak, after which the null gradually re- 
laxes towards its original potential configuration. Under 
continuous driving, it is unclear whether a true steady 
state would be set up, or whether the current sheet would 
continually grow in dimensions and intensity [see Ref . 77] . 
This is an open question for future study. For the case 
of transient driving, the peak current and reconnection 
rate increase linearly with the driving velocity. Here we 
examine more closely the sheet geometry, and its scaling 
with the driving velocity, and also investigate the scaling 
of this geometry, the peak current and peak reconnection 
rate with resistivity. 

As previously noted, the current sheet that forms in 
spine-fan reconnection is focussed at the null, locally 
spanning both the spine and fan. The sheet has a ten- 
dency to spread along the fan depending on the param- 
eters chosen (spreading is enhanced by lowering vq or 
increasing the plasma-/?). We examine four spatial mea- 
surements (see Figure 8b) associated with the current 
sheet, focussing on the time when the current magnitude 
is a maximum and defining the boundary of the sheet to 
be an isosurface at 50% of |j|ma£c- The sheet thickness is 
I, the length Ltot is the total extension in the yz-plane 
(normal to j), Lc is the length of the 'collapsed' section 
(within which the sheet contains both spine and fan), 
and the width w is the extension of the sheet along the 
x-direction (parallel to j). 

The scaling of these dimensions with (peak) driving ve- 
locity Vq is shown in Figure 9 (we fix 77 = 5 x 10~^). The 
angle 9 between the current sheet and the z — plane 
can be seen to increase as the driving velocity increases. 
This can be put down to the fact that the stronger driv- 
ing creates a stronger Lorentz force — the force that is 
responsible for the collapse. As expected, Lc increases as 
Vq increases. This is a result of the fact that the spine 
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footpoints are sheared further for larger vq, and there 
exists in fact a close correspondence; {LcCOs9)/2 ~ vqt. 
In contrast to Lc, Ltot shows a linear decrease with vq 
(as does w, see Ref. [77]), showing that as the collapse 
becomes stronger the distortion of the magnetic field fo- 
cusscs closer and closer around the null itself. The de- 
cline in Ltot with increasing vq must of course cease once 
Ltot = Lc, as is the case for the strongest driving consid- 
ered. Examining finally the sheet thickness /, any vari- 
ation is within the error bars of our measurements, and 
moreover the resolution is not sufficient for firm conclu- 
sions to be drawn. 

We turn now to consider the scaling of the current 
sheet with 77, setting vq = 0.02, see Figure 10. As 77 de- 
creases, jmax increases, while the reconnection rate de- 
creases. In both cases, with the limited data of this pre- 
liminary study, the proportionality appears to be some- 
where between power law and logarithmic. That the run 
with the largest resistivity does not seem to fit the trend 
for the reconnection rate is likely to be because the cur- 
rent significantly dissipates before reaching the null itself 
due to the high resistivity (77 = 0.002). Accompanying 
the increase in jmax with rj is, as expected, a decrease 
in the thickness /. On the other hand, the overall di- 
mensions of the sheet. Ltd and w, seem to be unaffected 
by rj, to within our measurement accuracy. Finally, as 
77 decreases and the current becomes more intense, the 
collapse becomes more pronounced as evidenced by in- 
creases in both Lc and 0. 

The relationships discussed briefly above certainly war- 
rant further investigation with carefully designed, higher 
resolution simulations, as do the corresponding scalings 



for the continuously driven case. 



VI. CONCLUSION 

We have here outlined a new categorisation of recon- 
nection regimes at a 3D null point. In place of the two 
previous types, namely, spine and fan reconnection, we 
suggest that three distinct generic modes of null point 
reconnection are likely to occur. The first two are caused 
by rotational motions, either of the fan or of the spine, 
leading to either torsional spine reconnection or torsional 
fan reconnection. These involve slippage of field lines in 
either the spine or the fan, which is quite different from 
classical 2D reconnection, but does involve a change of 
magnetic connection of plasma elements. 

Even though pure spine or fan reconnection may occur 
in special situations (such as when V ■ v = or there are 
high-order currents), it is much more likely in practice 
that a hybrid type of reconnection takes place that we 
refer to as spine-fan reconnection. This is the most com- 
mon form of reconnection that we expect to see in three 
dimensions at a null point, since it is a natural response 
to a shearing of the null point. It is most similar of all the 
3D reconnection regimes to classical 2D reconnection and 
involves transfer of magnetic flux across both the spine 
and the fan. It possesses a diffusion region in the form of 
a current sheet that is inclined to the fan and spine and 
has current localised in both the spine and fan, focussed 
at the null. 

In future, much remains to be determined about these 
new regimes of reconnection that have been observed in 
numerical experiments. One is the shape and dimensions 
of the diffusion regions and their relation to the driv- 
ing velocity and the magnetic diffusivity. Another key 
question is: what is the rate of reconnection at realis- 
tic plasma parameters, and is there a maximum value? 
Since the analytical theory is so hard in three dimen- 
sions, progress is likely to be inspired by future carefully 
designed numerical experiments. 
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APPENDIX A: ALTERNATIVE TORSIONAL 
FAN SOLUTION 



Here we present another pure nonideal solution for tor- 
sional fan reconnection. It possesses a much simpler elec- 
tric current, but the assumed form for the magnetic dif- 
fusivity has to be more complex and vanish on the fan. 
We follow section IV but consider a different form for the 
magnetic field of a double-spiral point, namely. 



{Br, B^, Bz) 



Bo 



R, 



jpzR 
2b ' 



-2z 



(Al) 



and corresponding current components 



Bojg 
HbLo 



(-ii?, 0, 



The field line equations for a mapping from an initial 
point (i?o, 7 b) to any other point (i?, 0, z) are 



and the inverse mapping is 



i?o = i?e-^, 0o-0-iio(l-e-2S), (A3) 



where S = — ilog(z/6). 

Assuming, as before, that z) vanishes in regions 
(1) and (2), we evaluate it in region (3) of Figure 6 by 
integrating from a point T{R, b) on the top of the diffu- 
sion region to a point P{R, z) inside the diffusion region. 
After using Equation (7) and the mapping (A2) and set- 
ting <I>e = 0, the expression for the potential at P[R, z) 
then becomes 



nidO 



Ab^rio 
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As an example, we adopt the following form for the 
magnetic diffusivity inside the diffusion region (D) 



which peaks above and below the null point and vanishes 
on the boundary of D. We have also chosen it to vanish 
on the fan plane since jz vanishes there and in order to 
facilitate a closed form solution with continuous physical 
quantities. After substituting into (A4) and using the 
mapping (A2) and inverse mapping (A3), we find the 



potential throughout the diffusion region as 



<^>md{R,z) = 
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4(m -2n + 2)6"'/2+3am 4(„j - 2n + 2)b"+^a' 
and the corresponding components of electric field as 
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We note that the solutions in the lower half plane z < 
may be obtained simply by replacing S = — ilog(z/6) 
by S' = — ^log(— z/6). There is a term in E^ that be- 
haves like R^/b^ and so is usually discontinuous at the fan 
plane. This discontinuity may, however, be avoided by 



15 



setting TO = 2 and balancing it with the term in z^^'^W^ 
when a^l\? = 12(n - l)(n + 6)/[5(n - 4)(n + 1)]. As a 
simple example, let us consider n = 6 and — 726^/7. 
Then, in order to find z) in region (4) of Figure 
6, as before, we calculate its value at any point Q{a,Zs) 
where the field line leaves the diffusion region, and then 
project that value along the field line. Thus, after putting 
(R,z) = {a,Zs) into Equation (A5), we obtain 



nidO 



1 Zg 


171.? 


48 ^ 76 


280fe2 


29zf 


z! 1 


^ 12066 


^ 1668 



As before ideal MHD holds in region (4). and so 



^nid{R, z) can be calculated from the fact that it is con- 
stant along the field line [zE? = ZgO^) joining Q to P. 
The electric field components follow and can be used to 
calculate the reconnection rate and the flow velocity in 
the usual way. 

Comparing the two forms of torsional fan solution, 
the advantage of the one presented in Section 4 is that 
the magnetic diffusivity peaks at the null point, but the 
slightly unwelcome feature is that the current density 
vanishes along the spine and fan. By comparison, the 
solution in this appendix has a current that vanishes at 
the null point but nowhere else in the spine or fan, but 
its disadvantage is that we had to choose the diffusivity 
to vanish in the fan. 
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